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ABSTRACT 


The  detection  capability  of  a  global  seismic  network  is  examined  on  the  basis  of  a 
probabilistic  model:  Given  the  location  of  se'smograph  stations  with  known  background 
noise  level,  a  worldwide  grid  of  epicenters,  and  amplitude-distance  attenuation  curves,  the 
detection  capability  is  expressed  by  the  magnitude  corresponding  to  a  fixed  probability 
that  a  specified  minimum  number  of  stations  detect  an  event. 


A  globally  distributed  network  of  50  stations  is  selected  and  multiwave  detection  cri¬ 
teria  are  applied  which  take  variation  of  attenuation  for  different  wave  types  into  con¬ 
sideration.  The  extension  of  amplitude-attenuation  curves  to  include  core  phases  is  inves¬ 
tigated  and  effects  of  regional  attenuation  are  studied.  Depending  on  ihe  detection  cri¬ 
terion  and  attenuation  curve,  magnitude  thresholds  of  a  50-station  network  can^ignifi- 
cantly(vary.  "w> , 


,r'ti 


*-'-’7  The  largest  influence,  however,  results  from  changes  in  station  noise  data.  Reliable 
estimates  of  noise  statistics  can  only  be  derived Trorn  continuous  measurements  in  connec¬ 
tion  with  routine  event  readings.  An  indirect  method  to  calculate  noise  parameters  is 
demonstrated  by  using  station  detection  thresholds.  After  estimation  of  noise  statistics 
from  station  reportings  to  the  International  Seismological  Centre  -f!SG)  for  1980, 
predicted  magnitude  detection  thresholds  and  empirical  values  are  in  good  agreement.  c; 


Magnitude  thresholds  for  teleseismic  detection  range  between  4.0^771^4.5  in  the 
northern  hemisphere  and  between  4.5^771^4.9  in  the  southern  hemisphere.  Model  calcu¬ 
lations  show  that  these  thresholds  can  be  lowered  at  least  half  a  magnitude  unit  after 
current  seismological  instrumentation  and  observatory  practice  have  been  improved. 


I. 


Introduction 


"Although  seismological  capacity  for  identifying  underground  nuclear  explo¬ 
sions  may  now  be  secondary  to  the  political  will  of  parties  engaged  in 
Comprehensive  Test  Ban  negotiations  it  is  still  important  to  present  the 
clearest  possible  evaluation  of  the  role  seismology  might  play  should  a 
Comprehensive  Test  Ban  become  reality." 

This  quotation  from  a  12  year  old  paper  (Marshall  and  Basham,  1972)  is  still  an 
adequate  description  of  the  general  purpose  of  studies  on  seismic  verification  of  nuclear 
test  ban  treaties.  In  a  more  specific  sense  we  want  to  assess  the  detection  threshold  of  a 
network  of  modern  seismic  stations.  Detection  thresholds  will  be  given  in  terms  of  magni¬ 
tude.  Therefore,  thresholds  described  herein  apply  to  both  shallow  earthquakes  and 
underground  explosions  without  regard  to  source  type.  The  important  questions  of  source 
identification  and  yield  estimation  are  not  addressed  in  this  paper.  Previous  detection 
studies  include  the  SIPRI-report  (Davies,  1968),  an  analysis  initiated  by  the  United 
Nations  (Basham  and  Witham,  1970),  and  a  report  of  a  group  of  seismologists  to  the 
Conference  on  Disarmament  in  Geneva  (CCD/558,  1978).  All  of  these  assessments 
present  conceptually  similar  schemes  whereby  worldwide  existing  seismological  facilities 
are  applied  to  a  straightforward  statistical  estimation  model:  Given  a  globally  distri¬ 
buted  network  of  seismograph  stations  with  known  background  noise  level,  a  set  of  epi¬ 
center  locations,  and  standard  amplitude-distance  attenuation  curves  for  seismic  waves, 
the  detection  capability  is  expressed  by  the  magnitude  corresponding  to  a  fixed  probabil¬ 
ity  that  a  specified  minimum  number  of  stations  detect  an  event. 

A  second  group  of  papers  (Kelly  and  Lacoss,  1969),  (Report  US/GSE/7,  1980) 
describe  a  different  approach  to  examine  the  detection  capability  of  a  seismic  network  by 
including  average  worldwide  seismicity.  Using  known  earthquake  recurrence  rates  a  syn¬ 
thetic  list  of  events  is  produced  as  a  reasonable  approximation  to  those  actually  observed 
in  a  specific  time  interval.  Keeping  station  parameters  unchanged,  this  method  allows  an 
independent  check  on  the  results  of  studies  of  the  first  kind. 

Finally,  a  third  procedure  starts  from  real  data  collected  during  special  experiments 
(Lacoss  et  al.,  1974)  or  published  in  bulletins  by  international  agencies  (e.g.,  International 
Seismological  Centre  (ISC)  in  Newbury,  UK)  and  estimates  detection  thresholds  with  the 
use  of  gaussian  or  maximum-likelihood  techniques  (Ringdal  et  ah,  1977),  (Ringdal,  1984). 

Detection  capabilities  of  seismic  networks  differ  substantially  as  a  result  of  these 
various  approaches,  in  general  the  detection  threshold  increases  in  the  order  of  the 
described  procedures.  Some  reasons  for  these  differences  are  obviously  due  to  the  differ¬ 
ence  between  operational  station  performance  used  in  the  last  approach  and  the  idealized 
assumptions  based  on  pure  noise  statistics  which  are  input  to  the  process  mentioned  first. 
Smaller  discrepances  simply  reflect  the  difficulty  in  making  this  type  of  estimate  and 
should  be  kept  in  mind  in  judging  on  the  accuracy  of  the  results. 


This  report  has  three  main  purposes: 


(1)  To  examine  recently  published  approaches  to  network  detection 
capability  estimation  by  using  multiwave  detection  criteria.  A  computer 
coded  version  of  this  procedure  (Ciervo  et  ah,  1983)  was  made  availab’e 
at  the  Center  for  Seismic  Studies  (CSS)  in  Arlington,  VA.  This  code, 
called  "Seismic  Network  Assessment  Program  for  Detection  (SNAPD)," 
not  only  models  the  propagation  of  p-waves  which  were  employed  in  pre¬ 
vious  programs  (Wirth,  1977)  but  takes  all  relevant  seismic  phases  into 
account  and  calculates  wave  attenuation  and  travel  time  as  a  function  of 
regional  media  characteristics  and  event  type. 

(ii)  To  study  the  influence  of  geophysical  input  parameters  on  the  out¬ 
come.  These  parameters  include  especially  the  amplitude  attenuation 
curves  at  teleseismic  distances  and  also  the  extension  to  core  phases. 
Variations  of  attenuation  in  tectonic  and  stable  areas  at  regional  dis¬ 
tances  are  of  special  importance  for  detailed  epicenter-staticr.  configura¬ 
tions.  Special  attention  will  be  given  to  amplitude  attenuation  curves 
derived  from  seismograms  of  underground  nuclear  explosions;  otherwise 
earthquake  data  arc  included  using  shallow  events. 

(iii)  To  reconsider  the  results  with  most  recent  station  data  which 
became  available  after  installation  of  a  large  number  of  digitally  recording 
seismographs.  Some  conclusions  can  already  been  drawn  from  the  data 
reported  to  the  ISC  for  1980.  These  data  are  used  to  compare 
maximum-likelihood  detection  estimates  (Ringdal,  1984)  with  results 
from  the  new  network  detection  code. 

II.  Analysis  of  Network  Detection  Probabilities 

The  statistical  model  and  basic  computational  procedures  are  described  in  this  sec¬ 
tion.  The  model  includes  various  parameters  to  be  known  at  the  beginning.  The  most 
important  are: 

•  Seismic  station  locations  and  their  noise  statistics  (mean  and  variance) 

•  Amplitude  distance  relations  for  several  phases  (besides  p-waves  and  R- 
waves  the  prominent  regional  phases  pg  and  lg  are  used) 

•  Signal-to-noise  ratio  required  for  detection. 

The  procedure  is  then  to  do  the  following: 

•  Select  a  source  location  and  compute  detection  probabilities  for  each  station 
as  a  function  of  event  magnitude. 

•  Find  the  lowest  magnitude  for  which  there  is  a  defined  probability  to  meet 
the  specific  detection  criterion.  (The  most  frequen*  criterion  for  global 
detection  studies  asks  for  a  90  per  cent  probability  of  detecting  p-waves  by 
at  least  four  stations.) 


First,  a  single  station  is  considered  and  the  probability  that  it  detects  a  certain  wave 
is  derived.  Then  the  multiwave  detection  probability  and  the  network  capability  are 
defined.  We  closely  follow  the  notation  by  Ciervo  et  ad.,  (1983).  Further  details  can  be 
found  in  Wirth,  (1977)  and  Elvers,  (1980). 


p.jl^  denotes  the  probability  that  wave  k  propagated  from  epicenter  j  will  be  observed  at 
station  i.  It  is  givin  by 
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where 


(1) 


(2) 


is  the  normal  cumulative  probability  function. 

In  Equation  (l)  signal  and  noise  are  assumed  to  be  log  normally  distributed  and  the 
log  of  the  noj.se  amplitude  has  expectation  Hn  and  variance  <7^  and  the  log  of  the  signal 
variance  is  o * .  <7^  defines  the  additional  variance  of  the  log  signal  amplitude  for  wave  k, 
given  an  m ^  value;  hence  (7^=0  if  k  denotes  the  p-wave.  A  station  i  is  supposed  to 
detect  wave  k  provided  that  the  ratio  of  signal  to  noise  is  at  least  r  .^. 

Given  an  event  at  epicenter  j  of  magnitude  rm  and  distance  A^.  from  station  i,  the 
amplitude  of  v/ave  k  at  that  station  is  calculated  as 

lo8-^  =  U) ( A,'y)  a)log(^,y)  ']c  (3) 


for  both  stable  (a=S)  and  tectonic  ( a=T )  media.  For  waves  other  than  conventionally 
used  p-phase  m ^  has  to  be  converted  into  an  adequate  magnitude  m ^  given  the  respective 
regression  formula 

mki=z\KE\k^K^kmb  (4) 


The  b ^  and  in  Equation  (3)  are  attenuation  table  entries,  and  is  the  epicenter- 
station  calibration  term  for  wave  k. 


If  wave  k  does  not  require  regional  attenuation  or  if  A>25 deg  then  logA^  is  com¬ 
puted  directly  from  (3)  using  a  stable  medium  attenuation  table.  Otherwise  for  regional 
distances  (AsC25de£f) 

log ^ ijk  ~  ( 1  ~  wij ) lu gA ijk  "*■  wi^°&A  (5) 

where  W-.  is  the  regional  path  weight;  i.e.,  the  ratio  of  the  length  of  the  wave  path  in  tec¬ 
tonic  meiia  to  the  total  great-circle  path  length  A t  Especially  if  the  epicentral  path  is 
assumed  to  have  passed  a  region  that  severely  attenuates  lg-waves,  then  logA;  <.  =—oc. 


If  the  attenuation  table  entries  are 


then  for  6,  . <A--<6 , 
it — 1  ij  r. 


6A/=JjfcandcA<HcJt 


except,  that  if  c^=0,  linear  interpolation  is  used  for  b: 


^rh-(VVth-h-ii/(Vh-.> 


The  station’s  probability  of  detection  given  by  Equation  (1)  is  influenced  by 

the  reliability  R.  of  its  operation,  p includes  therefore,  a  factor  R^,0<R^l.  R^  is  of 
course  dependent  on  a  number  of  local  circumstances  which  are  not  well  known.  Usually 
we  set  J?.=l,  assuming  perfect  operation. 

Given  the  probability  for  a  single  station  to  detect  an  individual  wave,  we  have  to 
develop  a  procedure  for  multiwave  network  detection  criteria  that  use  combinations  of 
dependent  wave  arrivals  at  individual  stations.  An  essential  feature  in  the  development  of 
the  model  is  the  assumption  that  a  minimum  of  four  phases  (not  more  than  two  of  which 
are  recorded  at  the  same  station)  are  required  for  detection  assessment.  A;  hough  relaxa¬ 
tions  of  this  requirement  are  possible,  they  do  not  seem  to  be  desirable  with  regard  to 
later  location  procedures. 

As  an  example  we  illustrate  the  detection  criterion 


(p.n'J/2  u  p/4 


which  literally  means  that  a  network  detection  consists  of  at  least  a  two-station  detection 
of  pg  and  lg  regional  waves  or  at  least  a  four-station  detection  of  p-waves.  Expression  (8) 
consists  of  two  subcriteria: 


D1  =  (p,  n  V  ft  D2  =  pji 


which  specify  different  wave  combinations. 

The  detection  criteria  such  as  Equation  (8)  have  to  be  decomposed  and  reduced  to  a 
set  of  canonical  probabilities.  Therefore,  the  logical  expression  (8)  is  transformed  into  an 
algebraic  expression  involving  the  marginal  probabilities  of  independent  subcriteria  and 
the  joint  probabilities  of  dependent  pairs  of  subcriteria.  This  transformation  eliminates 
all  logical  "or"  ((_j)  among  subcriteria  by  use  of  the  elementary  rule 

prob  [D1/n[jD2/m]  =  prob[D^/n  j  +  proi[D2/m]  (9) 

—prob  [D^npjD2/™] 

1  2 

where,  if  D  and  D  have  no  waves  in  common 

prob  [D1/r?p|D2/m  ]  =prob  [D1/n  ]  prob  [D2/rn  ]  (10) 


In  case  of  Equation  (10)  independent  calculations  for  each  individual  wave  define  the 
probability  proi>[D^-  ]  that  exactly  n  out  of  N  stations  detect  wave  combination 
D  .  If,  however,  the  detection  subcriteria  I)1  and  D  have  waves  in  common  like 


Wnpnup/* 


then  computation  of  the  joint  probability  (10)  does  not  split  into  independent  probabili¬ 
ties. 

Clearly,  the  probability  pro£>[D  / n]  of  at  least  n  detecting  stations  follows  as  : 


i  v 

p:ob  [D%  1  =■  E  Frob  [D^  ] 


Thus,  we  need  only  to  compute  the  probability  of  exactly  n  detecting  stations. 

In  a  concluding  step  a  binary  search  is  used  to  find  the  magnitude  value  that  results 
in  pt,  the  threshold  probability  for  network  detection.  Let  m .  be  a  sequence  ol  test  mag¬ 
nitudes  ^ mm5jf  such  that  ea°h  m-  results  in  a  network  detection  probability 

P--  Initially  If  P^Pt,  the  search  is  termimited;  otherwise  m2~mmax  and  if 

p2^pt,  the  search  is  also  terminated.  Assuming  p^<pt<p 2 
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)/2  *  — 1,2.3,.. 


where,  if  p-  ,  is  less  than  p ^  m  is  the  last  test  magnitude  for  which  the  corresponding 
network  probability  is  greater  than  p(,  and  vice  versa. 

* 

The  search  is  continued  until  »  =  »  is  reached  such  that 

I  P,- — Pfl 


corresponding  to  a  threshold  magnitude  m.- 


III.  Detection  Capability  of  a  Global  Network 

In  selecting  a  network  for  detection  of  seismic  events  baied  upon  existing  seismo¬ 
graph  stations,  it  is  desirable  to 

•  arrive  at  a  relatively  uniform  geographical  distribution  of  stations 

•  select  stations  with  modern  instrumentation  and  optimum  detection  capa¬ 
bilities. 

With  these  criteria  in  mind  a  selection  of  relatively  few  stations  is  considerably  more 
effective  than  using  all  (about  1000)  stations  that  routinely  report  to  one  of  the  interna¬ 
tional  data  centers.  In  a  previous  report  (CCD/553,  1978)  a  network  of  50  stations  has 
been  composed  which  was  judged  to  produce  the  best  seismological  results  currently 
achievable  for  teleseismic  dctectior  of  seismic  events. 

These  stations  whose  geographical  distribution  is  shown  in  Figure  1  and  whose  coor¬ 
dinate0  and  further  parameters  are  given  in  Table  1  arc  used  as  a  reference  network  in  our 
study.  The  global  station  distribution  is  not  as  uniform  as  desirable  having  36  stations  in 
the  northern  hemisphere  compared  to  14  stations  in  the  southern  hemisphere  or  31  sta¬ 
tions  in  the  eastern  hemisphere  compared  to  19  in  the  western  hemisphere,  but  it  reflects 
to  some  extent  the  distribution  of  land  masses  on  earth. 

Mean  value  and  standard  deviation  of  logarithmic  station  noise  which  arc  used  to 
calculate  station  detection  probabilities  from  Equation  (1)  are  also  included  in  Table  1 
(columns  4  and  5).  Due  to  the  quoted  report  (CCD/558,  1978)  average  noise  levels  were 
partly  derived  from  published  noise  power  spectra,  partly  estimated  from  magnification 
curves  of  seismographs.  Because  of  lack  of  adequate  measurements  an  adhoc  procedure 
was  used  to  estimate  the  variance  of  the  station  roise;  those  stations  with  higher  noise 
level  were  abo  assigned  greater  variance. 
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As  can  be  seen  from  Equation  (l)  it  is  the  combined  effect  of  variance  in  noise 
amplitude  and  signal  amplitude  which  influences  the  station  detection  threshold;  it 
decreases  as  the  denominator  in  Equation  (l)  is  increased  provided  that  the  detection  pro¬ 
bability  is  less  than  0.5,  and  vice  versa.  For  the  network  the  detection  threshold  generally 
decreases  v/hen  signal  or  noise  variance  is  increased.  For  the  50  station  network  -  intro¬ 
duced  in  Figure  1  and  Table  J.  -  a  constant  value  for  the  standard  deviation  of  signal 
amplitude  (0.2  in  logarithmic  units)  was  used.  Test  runs  showed  that  doubling  this 
parameter  to  0.4  result  in  a  very  small  dilference  of  the  network  detection  threshold  (not 
exceeding  0.1  magnitude  unit).  Consequently,  network  capability  is  not  very  sensitive  to 
this  parameter. 

There  are  two  other  input  parameters  to  Equation  (1)  which  are  to  be  assumed  in  an 
adhoc  manner:  the  reliability  factor  R  describing  grossly  the  station  operation  (up-time) 
is  set  in  our  calculations  to  1.0.  Earlier  studies  (Ringdal  et  al.,  1977)  have  estimated  this 
parameter  to  range  from  0.8  to  1.0  for  most  of  the  stations  selected  for  our  network. 

A  minimum  signal-to-noise  ratio  has  t.o  he  chosen  to  detect  seismic  signals  emerging 
from  background  noise.  Throughout  this  study  a  s/n  ratio  (  r.,  in  Equation  (1)  )  of  1.5 
was  chosen.  The  capability  results  are  easily  transformed  to  correspond  to  other  s/n 
ratios,  since  this  parameter  occurs  as  a  difference  (m  —  logr^)  in  Equ  ion  (1).  Thus,  a 
simple  relationship  exists  between  the  chosen  value  of  and  the  corresponding  magni¬ 
tude  level.  If  r.^  is  increased  from  1.5  to  3.0,  for  example,  the  threshold  will  be  increased 
by  [logZ  —  logl.5)“0.3  magnitude  units. 

a.  Amplitude-Attenuation  Curves  at  Teleseismic.  Distances: 

25dcg<  A<100dcg 

Several  investigations  of  amplitude  attenuation  with  distance  have 
been  made  since  the  pioneering  work  of  Gutenberg  and  Richter,  (1956). 
According  to  their  paper  we  will  summarize  amplitude-distance  relations 
from  Equation  (3)  in  the  form 

D  (A)  -  b  (A)  +  c -log  (A) 

So  we  can  interpret  B-values  in  terms  of  magnitude  units. 

In  Figure  2  the  results  of  Evernden  and  Clark,  (1970),  Veith  and 
Clawson,  (1972)  and  NORSAR  (  )  are  compared  to  the  Gutenberg  and 
Richter  curve.  Whereas,  Evernden  and  Clark  as  well  as  Veith  and  Clawson 
use  mostly  explosions  and  a  station  network  of  mainly  LRSM-stations  in 
the  U.S.,  the  NOIlSAR-curve  is  derived  from  observations  of  teleseismic 
events  at  a  single  array  site. 

All  curves  have  been  arbitrarily  connected  at.  regional  distances  tc 
focus  upon  the  differences  in  the  teleseismic  window.  The  principal  differ¬ 
ence  between  Veith-Clawson  and  NORSAR  curves  on  one  hand  and 


Gutenberg-Richter  and  Evernden-Clark  curves  on  the  other  hand  appears  in 
the  fact  that  the  former  smooth  amplitude  variations  from  mantle  discon¬ 
tinuities,  whereas  the  latter  indicate  several  step-wise  changes  in  amplitude 
as  a  function  of  distance.  The  difference  between  Veith-Clawson  and  NOR¬ 
SAR  attenuation  curves  is  nearly  a  constant  0.1  magnitude  unit  over  the 
whole  teleseismic  range.  This  may  result  from  the  fact  that  the  Veith- 
Clawson  curve  was  corrected  for  surface  focus  events;  in  comparison  the 
NORSAR  curve  comprises  measurements  from  shallow  earthquakes.  The 
obvious  distinction  between  the  curve  of  Gntenherg-Richter  and  Evernden- 
Clark  appears  in  the  far  teleseismic  portion  where  a  difference  of  0.3  magni¬ 
tude  units  can  be  found.  This  is  a  consequence  of  the  way  Evernden  and 
Clark  have  chosen  to  normalize  their  data. 

Some  of  these  differences  can  be  seen  in  the  network  detection  capabil¬ 
ity  which  is  shown  in  Figures  3-6.  Using  a  15  deg  epicenter  grid  the  90  per 
cant  probability  of  at  least  four  detecting  stations  was  calculated.  For  the 
NORSAR  attenuation  function  (Figure  3),  which  will  be  used  as  a  reference 
curve  in  this  study,  the  magnitude  threshold  is  estimated  to  be  from  3.4 
to  3.7  in  Europe  and  Scandinavia,  3.7  to  3.8  in  North  America,  Asia, 
and  Arctica,  rn^  3.7  to  4.0  in  Africa  and  most  parts  of  South  America, 
whereas  we  get  values  up  to  rn,  =4.5  in  the  Pacific  region.  The  slight  differ¬ 
ence  between  the  eastern  and  western  hemisphere  (about  0.2  magnitude 
units)  as  well  as  the  large  difference  of  more  than  one  magnitude  unit 
between  the  northern  and  southern  hemisphere  mainly  result  from  the  sta¬ 
tion  distribution  of  the  network.  The  high  station  noise  at  the  sites  in  the 
Pacific  (New  Zealand  and  Samoa)  gives  n  additional  contribution  to  the 
low  detection  capability  in  the  southern  hemisphere. 

As  expected  from  the  preceding  discussion  of  attenuation  curves 
(Figure  2)  detection  thresholds  for  the  network  increase  globally  by  0.1 
magnitude  unit  using  the  Veith-Clawson  data.  The  Evernden-Clark  curve 
has  a  remarkable  effect  in  lowering  the  detection  threshold  in  the  southern 
hemisphere  by  about  0.3  magnitude  units  (Figure  5),  again  not  unexpected 
from  the  shape  of  the  attenuation  curve  which  shows  pronounced  lower  13- 
values  at  far  teleseismic  distances  than  any  other  attenuation  cur.*;  in  Fig¬ 
ure  2.  Finally,  Gutenberg-Richter  curve  yields  detection  thresholds  (Figure 
6)  very  similar  to  Veith-Clawson  curve  (Figure  4). 

We  can  conclude  from  these  calculations  that  step-wise  changes  in  the 
attenuation  curve  do  not  significantly  effect  the  global  detection  capability 
of  a  50-station  network.  Of  course  these  discontinuities  have  remarkable 
focussing-defocussing  effects  for  specific  epicenter-station  configurations,  but 
these  arc  smoothed  and  can  not  be.  resolved  by  global  grids  of  15  deg  size. 
These  peculiarities  are  better  implemented  by  use  of  epicenter-station  cali¬ 
bration  factors  expressed  by  in  Equation  (3). 

For  all  computations  we  kept  the  probability  level  at  90  per  cent.  In 
changing  this  parameter  one  can  significantly  influence  detection  thresholds 


(of  course  the  seisriological  capability  i3  not  changed  at  all).  At  the  30  per 
cent  probability  level  (Figure  7)  we  get  0.2  -  0.4  magnitude  units  lower 
thresholds  compared  to  the  comnmnly  accepted  90  per  cent  probability 
level. 

Another  way  of  demonstrating  this  difference  is  to  calculate  the  net¬ 
work  detection  probability  for  a  fixed  magnitude  value.  Figure  8  shows  the 
global  distribution  of  probabilities  to  detect  a  magnitude  4  event.  Besides 
Antarctica,  New  Zealand,  the  Pacific  islands,  and  the  tip  of  South  America 
the  chance  to  detect  events  on  land  down  to  this  size  at  (at  least)  four  sta¬ 
tions  is  higher  than  80  per  cent. 

b.  Extension  of  Amplitude-Attenuation  Curves  beyond  100  rie; 
Distance 

The  most  prominent  result  of  the  last  section  is  the  clear  difference  in 
detection  capability  between  northern  and  southern  hemisphere  regardless 
what  type  of  attenuation  function  was  applied  in  the  teleseismic  window. 
To  reduce  the  detection  threshold  in  the  southern  hemisphere  a  substantial 
increase  of  the  number  of  stations  in  that  region  has  often  been  recom¬ 
mended.  Although  this  is  a  solution  in  principle,  there  are  several  practical 
problems  with  its  realization.  The  major  part  of  the  southern  hemisphere  is 
covered  by  deep  ocean  areas  where  installation  and  maintenance  of  seismo¬ 
graphs  is  still  difficult  and  expensive.  Islands  are  known  for  a  high  micrcse- 
ismic  noise  level  and  poor  detection  capability. 

An  attractive  alternative  for  improving  the  detection  capability  in  the 
southern  hemisphere  is  provided  by  signals  which  have  travelled  through  the 
earth’s  core  and  are  routinely  detected  at  stations  beyond  100  deg  distance. 
In  a  specific  window,  i.e.,  142  deg<&<  152  deg  these  refracted  core  phases 
offer  even  better  detection  possibilit'es  than  earlier  described  direct  p-waves. 
Numerous  studies  have  shown  that  observation  of  high  amplitudes  for  vari¬ 
ous  pkp-branchc3  can  be  a  powenul  tool  in  lowering  detection  thresholds. 
(Blandford  and  Sweetser,  1973),  (Quemar,  1973.)  This  can  easily  be 
demonstrated  with  a  seismogram  of  a  F  rench  .  dear  underground  explosion 
exploded  at  Mururoa  atoll  (21  S,140  W  and  recorded  with  the  Graefenberg 
array  in  Germany  (A ~  143.7 deg).  Figure  9  shows  a  recording  at  all  13 
vertical  elements  of  the  array  (for  a  more  detailed  description  see  Harjes 
and  Seidl,  1978)  and  the  beam-trace  on  top  from  which  a  displacement 
amplitude  of  4  nm  at  a  period  of  0.8  sec  was  measured. 

In  comparison  Figure  10  shows  recordings  of  the  same  event  by  the 
RSTN-stations  in  North  America.  These  stations  (Engdrhl  et  al.,  1982)  are 
new  borehole  installations  and  include  seismometers  (Geotech  S-750)  with 
high  sensitivity  in  the  short-period  band.  Traces  are  aligned  to  the  theoreti¬ 
cal  arrival  time  of  the  p-wave  which  is  marked  by  the  cursor  line.  The  sta¬ 
tions  have  a  distance  of  73  deg  to  88  deg  from  the  event.  Neither  the  origi¬ 
nal  (upper  part)  nor  the  narrow-band  filtered  (lower  part)  traces  meet  the 


detection  requirements  set  in  our  calculations.  Correspondingly,  this  event 
was  not  reported  by  inter  national  data  centers  which  restrict  their  jvent- 
defining  association  process  to  p-arrivals  within  100  deg  distance. 

Figure  11  shows  the  summary  of  station  reportings  for  this  event 
available  from  GTS/WMO-channels  at  the  CSS.  The  association  program 
implemented  at  the  CSS  found  the  questionable  event  by  using  pkp- 
observations  from  three  stations  in  Europe  (including  GRF).  With  the 
detection  criterion  we  applied  in  the  preceding  paragraph  (at  least  4  p- 
detections)  this  event  would  have  been  missed. 

Consequently,  we  amended  amplitude-distance  curves  beyond  100  deg 
as  shown  in  Figure  12.  There  are  two  curves  from  different  sources  (Bland- 
ford  and  Sweetser,  1973),  (Ringdal,  1984)  which  show  a  great  similarity 
although  Blandford  and  Sweetser’s  curve  is  based  on  a  much  more  general 
.ata  base  (ISC)  than  Ringdal’s  data  which  are  derived  only  from  NORSAR 
observations.  Detection  threshold  estimates  for  the  50-station  network  cal¬ 
culated  with  thes’  attenuation  curves  differ  less  than  O.i  magnitude  unit. 
Israelsson,  (1984)  recommended  the  requirement  of  at  least  one  p  detect’on 
additional  to  pkp-arrivals  to  define  an  event  to  avoid  large  location  errors 
because  the  stations  observing  pkp  might  be  clustered  in  a  narrow  geograph¬ 
ical  area.  Figure  13  gives  the  network  detection  capability  for  detecting  at 
least  4  p-  or  pkp-arrivals,  one  of  which  has  to  be  a  direct  p-wave,  i.e., 
obser\  at  a  station  within  less  than  100  deg  distance  of  the  event.  This 
result  c.  directly  compared  with  Figure  3  because  the  attenuation  func¬ 
tions  are  identical  up  to  100  deg  distance. 

The  inclusion  of  pkp-phases  yields  a  large  decrease  of  the  detection 
threshold  for  the  southern  hemisphere  (0.3  to  0.5  magnitude  units)  and 
divides  the  difference  in  detection  capability  between  northern  and  southern 
hemisphere  by  half  compered  to  the  only  use  of  p-arrivals  within  100  deg 
distance. 

Variation  of  Amplitude-Attenuation  Curves  at  Regional  Distances 

In  discussing  the  influence  of  varying  teleseismic  attenuation  functions 
on  the  network  detection  capability  we  always  used  the  same  curve  for 
regional  distances  (A<25  deg).  Although  the  teleseismic  part  (including 
pkp)  will  be  most  important  for  a  global  network  regional  differences  will 
have  some  effect  in  areas  (Europe  and  Scandinavia)  where  the  50-station 
network  is  highly  represented.  These  areas  give  us  the  opportunity  to  study 
the  effect  of  different  regional  attenuation  functions  which,  in  a  gross  sense, 
represent  ''stable"  (i.e.,  high  Q)  and  "tectonic"  (i.e.,  low  Q)  provinces. 

The  distinction  was  initially  introduced  for  the  North  American  con¬ 
tinent  taking  into  account  ths  different  crustal  structure  in  western  and 
eastern  U.S.  (Evernden,  1967).  In  this  section  we  discuss  only  the  influence 
of  variations  of  p  -  wave  which  is  seen  as  the  first  arrival  on  regional 


seismograms  recorded  at  distances  greater  than  1  deg.  To  emphasize  the 
difference  rather  extreme  representatives  of  published  attenuation  curves 
have  been  used,  namely  Evernden’s  "8.5"  curve  (Evernden,  1957)  derived 
from  data  in  the  eastern  U.S.  (b=-0.83,  c=-2  in  Equation  (3))  compared  to 
a  attenuation  curve  d  .  1  from  data  in  southwestern  U.S.  (Der  et  al., 
1982)  which  yield  b=0.1w-.  and  c=-3.803.  Thus,  tl  main  difference 
between  these  two  types  of  regional  attenuation  curves  is  that  the  "stable" 
p  drops  off  as  the  square  of  the  distance  while  the  "tectonic"  curve 
decreases  more  rapidly,  almost  with  the  fourth  power  of  the  distance. 

Another  peculiarity  at  regional  distances  is  thi  relative  maximum  in 
the  amplitude-distance  curve  as  a  result  of  the  20  deg  discontinuity  which  is 
observed  worldwide  with  differing  prominence  (Gutenberg  and  Richter, 
1956),  (Veith  and  Clawson,  1972).  To  emphasize  also  the  effect  of  the  20 
deg  discontinuity  the  Veith-Clawson  curve  has  been  appended  to  Evernden’s 
curve  in  the  distan  t  range  17  deg< A<25  deg.  On  the  opposite  side  the 
"tectonic"  •'.ttenuation  c\  rve  has  been  smoothly  connected  to  teleseismic 
distances  without  consideration  of  a  20  deg  discontinuity. 

Our  final  regional  extensions  of  attenuation  functions  are  plotted  in 
Figure  14.  It  is  again  noted  that  these  curves  are  artificial  compositions  to 
show  the  most  pronounced  effect  on  network  deter tion  capability.  The 
result  can  be  seen  in  Figure  15  and  Figure  16  which  are  to  be  compared  to 
Figure  13.  Using  specified  regional  attenuation  curves  generally  increases 
the  influence  of  station  distribution.  Global  differences  of  detection  thres¬ 
holds  are  pronounced  in  Figure  15  which  shows  the  influence  of  a  "stable" 
(Evernden  "8.5"+  20  deg  disc.)  regional  attenuation  curve.  At  close-to- 
station  distances  detection  thresholds  are  decreased  by  up  to  0.2  magnitude 
units  as  compared  to  Figure  13.  On  the  contrary  detection  thresholds  are 
increased  for  most  parts  of  Eurasia  and  America  by  up  to  0.2  magnitude 
units  for  a  "tectonic'1'  attenuation  curve  (Figure  16).  Worldwide,  this 
attenuation  curve  smoothes  the  difference  in  detection  capability  between 
northern  and  southern  hemisphere,  certainly  a  result  of  our  station  distribu¬ 
tion. 


In  summary  regional  variation  of  amplitude  attenuation  curves  can 
change  the  detection  threshold  of  a  global  50-station  network  by  as  much  as 
±0.2  magnitude  units. 


IV.  Regional  Detection  Probabilities  of  a  Global  Network 

Although  this  paper  mainly  deals  with  teleseismic  detection  capabilities  it  is  interest¬ 
ing  to  investigate  to  what  extent  the  magnitude  threshold  of  a  global  network  will  be 
influenced  by  including  other  than  p^  phases.  It  should  be  mentioned  at  the  beginning 
that  this  section  is  intended  as  an  amendment  to  teleseismic  capabilities.  To  study  the 
full  potential  of  regional  waves  for  detection  purposes  a  regional  station  network  has  to  be 
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introduced.  This  is  beyond  the  scope  of  this  report  and  for  details  on  this  matter  we  refer 
to  excellent  recent  review  articles  (Pomeroy  et  ah,  1982),  (Blandfi  ;d,  1981). 
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It  is  well  known  that  the  largest  amplitudes  on  a  Sort-period  regional  seismogram 
occur  within  the  lg-wavetrain  which  may  be  interpreted  as  a  superposition  of  a  large 
number  of  higher  mode  Love  and  Rayleigh  waves  -  lg-amplitudes  can  be  ten  times  larger 
than  the  maximum  amplitude  of  the  first  arrival  (p-wave)  at  the  same  distance  and  in  the 
same  period  band  around  one  sec.  The  actual  amplitude  can  drastically  vary  due  to  local 
geology.  An  extensive  literature  exists  on  amplitude-distance  attenuation  cu’-ves  for  lg  for 
various  regions  of  the  world  (for  reference  see  the  above  mentioned  review  articles). 

Again,  we  extract  two  extreme  representatives  of  published  curves  to  examine  the 
effect  on  ths  detection  capability.  As  representing  "stable"  regions  with  low  attenuation 
we  choose  an  attenuation  curve  published^VgNuttli,  (1973).  Nuttli  derived  an  amplitude 
decay  with  distance  proportional  to  A  '  (0.5deg< A<40deg)  corresponding  to  the 
shape  of  well-known  "Prague"-formula  (Vanek  et  ah,  1962)  adopted  by  LASPEI  to  be 
used  for  teleseismic  Rayleigh  wave  observations.  However,  Nuttli  derived  his  curve  from 
observations  of  onel  sec  lg-waves  in  the  eastern  U.S.,  whereas  the  "Prague"-formula  is 
valid  for  Raleigh  waves  around  20  sec. 

For  crustal  structures  in  "tectonic" ^provinces  like  the  western  U.S.  an  amplitude- 
distance  decay  for  lg  proportional  to  A  was  observed  (Der  et  ah,  1982).  The  selected 
attenuation  curves  are  plotted  in  Figure  17.  In  estimating  detection  capabilities  by  use  of 
lg-waves  it  has  to  be  mentioned  that  its  amplitude  is  not  only  regionally  varying  but  that 
it  also  can  be  totally  squelched.  Representing  a  wave  guide  phenomenon  lg-propagation  is 
seriously  effected  by  variations  in  crustal  thickness  (mountains, ocean-continent  and/or 
tectonic  boundaries).  To  take  these  inefficient  wave  paths  into  account  a  much  finer  than 
the  15  deg  x  15  deg  grid  used  in  our  calculations  has  to  be  applied. 

We  want  to  examine  whether  inclusion  of  lg-waves  has  an  effect  on  a  global  network 
at  all.  If  there  appear  regional  improvements  of  the  detection  capability  these  have  to  be 
verified  by  considering  the  corresponding  regional  crustal  structure.  Because  we  want  to 
insist  on  teleseismic  detections  we  use  lg-detections  only  if  at  least  one  station  of  the  net¬ 
work  observed  a  teleseismic  p-wave  (  A>25  deg).  So  we  required  at  least  two  lg- 
detections  and  1  teleseismic  p-detection  or  four  p-detections.  As  lg  is  observed  on  all 
three  components  of  ground  motion  an  azimuth  estimate  can  be  calculated  and  two  sta¬ 
tions  t»ro  sufficient  to  roughly  associate  the  event  origin.  An  event  is  declared  if  at  least 
one  p-observation  at  teleseismic  distam-s  confirms  this  association. 

Figure  18  and  Figure  19  show  the  detection  capability  using  this  criterion  for 
Nuttli’s  and  Der  et  al.’s  attenuation  function  respectively,  in  comparison  to  Figure  15  and 
Figure  16  which  show  the  corresponding  p-wave  detection  results  thresholds  are  lowered 
by  0.1  to  0.3  magnitude  units.  Again,  it  should  be  emphasized  that  these  improvements 
are  irrelevant  if  they  occur  in  oceanic  areas  because  lg-waves  disappear  after  crossing 
approximately  100  km  of  oceanic  structure.  Restricting  the  evaluation  to  continental 
areas  only  it  can  be  seen  from  Figure  18  and  Figure  19  that  the  detection  threshold  is 
mostly  influenced  in  the  southern  hemisphere  where  occasionally  two  stations  are  located 
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at  regional  distances.  In  Europe  and  Scandinavia  we  reach  already  a  high  capability  for 
four  p-detections  which  is  not  significantly  improved  by  the  additional  (two  lg  and  one 
telescismic)-detection  probability.  Finally,  it  might  be  mentioned  that  the  detection 
difference  caused  by  the  difference  of  the  two  attenuation  curves  (Figure  17)  is  only  mar¬ 
ginal. 

V.  Station  Data 

After  discussing  the  influence  of  amplitude-distance  curves  on  network  detection 
capability  the  remaining  factor  which  determines  detection  thresholds  is  the  noise  statis¬ 
tics  of  stations.  Published  magnification  curves  or  spectral  estimates  from  short  time 
windows  formed  the  basis  to  estimate  these  values.  More  reliable  data  can  be  obtained  by 
actually  measuring  the  performance  of  each  station  over  an  extended  period  of  time. 
Unfortunately  noise  measurements  are  not  reported  regularly  by  stations,  so  an  indirect 
method  was  applied  to  estimate  those  from  station  detection  thresholds.  The  data  base  of 
the  International  Seismologica!  Centre  (ISC)  for  1980  was  used  as  a  reference  system  and 
station  detection  values  computed  by  Ringdal,  (1984)  were  taken  to  derive  station  noise 
values.  The  following  procedure  was  used: 

•  given  :  station  detection  threshold  log  A/T 

•  calculate  :  amplitude  of  detection  threshold  A  at  T  =  one  sec 

•  assume  :  signal-to-noise  ratio  r  for  detection 

•  deduce  :  station  noise  amplitude 

Before  discussing  the  result  of  this  procedure  it  has  to  be  verified  that  the  selected 
stations  reported  to  the  ISC  in  1980.  In  column  6  of  Table  1  the  number  of  teleseismic 
detections  published  by  ISC  is  listed  for  all  stations  belonging  to  the  50-station  network. 
Of  course  no  noise  estimate  can  be  deduced  for  stations  which  do  not  appear  in  this  list. 
Fortunately,  for  some  stations  we  could  find  a  nearby  substitute.  They  can  be  recognized 
from  Table  1  having  additional  station  names  in  brackets  which  were  the  original  ones 
selected  in  our  50-station  network.  Finally,  we  were  able  to  calculate  "operational"  noise 
statistics  (MLE-SP  noise  in  Table  1)  for  40  stations  which  can  be  seen  to  differ  substan¬ 
tially  from  earlier  estimated  figures  (column  4  and  5  of  Table  1).  The  difference  can 
easily  be  demonstrated  in  Figure  20  which  shows  a  histogram  of  noise  data  we  used  in 

network  detection  studies  as  thick  lines  compared  to  our  new  values  (dotted  lines.) 

derived  from  Ringdal’s  maximum-likelihood  estimates  (MLE)  of  station  detection  thres¬ 
holds.  The  "operational"  noise  values  represent  much  closer  a  normal  distribution  peaked 
around  0. 8-0.9  (corresponding  to  6  nm  -  8  nm)  than  the  original  data  which  showed  12 
stations  with  a  mean  noise  value  of  2  nm.  It  might  be  mentioned  that  the  shilt  of  very 
sensitive  stations  (like  LAO  and  NAO)  is  biased  by  taking  the  ISC-bulletin  as  reference 
because  part  of  their  detections  do  not  meet  the  event  defirf  criteria  of  the  center.  But 
this  restriction  of  our  procedure  has  little  influence  on  the  owing  discussion  as  the  net¬ 
work  capability  is  largely  determined  by  the  general  increase  of  noise  values. 

To  elaborate  further  on  this  difference  of  noise  estimates  we  examined  data  from 
oui  German  station  at  Graefcnberg  (GRF).  We  could  use  an  adequate  data  set  while 


GRF  participated  in  a  Common  Data  Base  Experiment  (CDBE)  which  was  conducted  by 
Sweden  under  the  auspices  of  the  Geneva  group  of  seismologists  during  the  time  period 
Oct  1  -  Oct  15  1980  (Barkeby  et  al.,  1981).  Figure  21  shows  a  power  spectral  density 
(PSD)  plot  of  noise  measured  at  GRF  during  the  time  period  of  that  experiment.  This 
spectrum  represents  an  average  over  20  noise  samples  taken  in  time  intervals  of  four  min 
between  00:00  h  and  01:30  h  on  Oct  10.  Each  noise  sample  was  computed  from  a^ime 
series  of  2048  values  (At=0.05sec)  multiplied  uy  a  cosine-window.  We  read  0.4nm  j Hz 
at  1  Hz  from  this  PSD-plot.  Assuming  normal  distributed  noise  we  calculate  an  average 
zero-to-peak  noise  amplitude  using  following  formula  (Taylor,  1981). 

averageO-p-amplitude  ~  [2 BW(Hz)-$(/)(a'2/Hz)] 


where  BW  is  the  bandwidth  of  the  recording  system  and  $^/)  is  the  power  spectral  den¬ 
sity.  With  BW  =  5  Hz  for  GRF-seismograph  and  3>=0.4nm  I  Hz  we  get 

averageO—p— amplitude  ~  2 nm 


a  value  which  in  fact  was  supplied  for  network  capability  calculations  in  1978  (CCD/558) 
as  can  be  seen  from  Table  1  (column  4). 

During  the  Common  Data  Base  Experiment  stations  were  asked  to  report  noise 
measurements  in  addition  to  standard  parameters  and  were  instructed  to  read  the  max¬ 
imum  noise  value  within  30  sec  preceding  the  p-wave  (in  the  frequency  band  of  the  sig¬ 
nal). 


Figure  22  shows  a  histogram  of  these  noise  amplitudes  (ITanka  and  Henger,  1980)  for 
station  GRF.  Although  this  data  set  is  rather  limited  (70  values),  and  not  well 
represented  by  a  normal  distribution,  a  calculation  of  mean  ^nd  standard  deviation  shows 
good  agreement  with  Ringdal’s  figure  (Table  1).  Besides  i  can  be  seen  from  Figure  22 
that  there  is  a  smaller  peak  of  measurements  at  about  0.3  (ihe  value  derived  from  PSD- 
data).  The  majority  of  data,  however,  is  shifted  about  0.5  to  higher  noise  values. 

If  we  assume  that  the  procedure  demonstrated  for  GRF  is  generally  applicable  then 
the  conclusion  can  be  drawn  that  average  0-p-amplitude  noise  values  (from  PSD- 
estimates)  are  not  appropriate  as  input  data  in  the  detection  calculations  but  rather  that 
the  actually  measured  maximum  0-p-values,  taken  over  many  samples,  are  required.  Of 
course,  this  conclusion  is  closely  related  to  the  specified  signal-to-noise  ratio  r  =  1.5  at 
which  stations  are  assumed  to  detect  an  event.  In  the  algoritl  ^Equation  l)  one  could 
always  compensate  for  systematically  low  noise  estimates  by  in  asing  s/n-ratio  to  r  =  5 
(for  example).  Commonly  a  value  of  r  =  1.5  is  more  acceptable  if  correct  noise  data  can 
be  found. 

Thirty-nine  stations  (marked  by  an  asterisk  in  Table  l)  were  used  to  compare  our 
results  from  the  "model  network"  with  real  station  data.  The  geographical  distribution  of 
these  39  stations  can  be  seen  in  Figure  23.  To  some  extent  this  reduced  network 


represents  a  mce  even  global  coverage  than  the  50-station  network.  At  first  we  exam¬ 
ined  the  effect  of  the  reduction  in  station  number.  Figure  24  shows  the  detection  capabil¬ 
ity  of  this  reduced  network  whicn  i3  directly  comparable  to  Figure  3  because  the  same 
attenuation  curve  and  same  station  data  were  used.  Although  we  have  earlier  reported 
from  previous  studies  that  a  smaller  network  of  high-quality  stations  is  preferable  to  the 
full  use  of  all  worldwide  existing  stations  the  small  effect  of  reducing  our  50-station  net¬ 
work  by  more  than  20  per  cent  is  still  surprising.  The  difference  of  detection  threshold  is 
at  most  0.1  magnitude  unit  and  therefore  negligible  within  the  numencal  precision  of  our 
results. 

A  significant  change  is,  however,  accomplished  if  the  new  station  data  (noise  statis¬ 
tics)  are  introduced  (Figure  25).  Magnitude  thresholds  are  increased  by  about  half  a 
magnitude  unit.  Ringdal  (personal  communication,  1984)  used  his  station  detection 
thresholds  and  recurrence  estimates  from  observed  earthquakes  in  1980  to  obtain  detecta¬ 
bility  thresholds  for  regions  with  sufficient  number  of  observations  (Kelly  and  Lacoss, 
1969).  His  results  for  the  same  network  are  shown  in  brackets  on  Figure  25  and  a  very 
good  agreement  to  our  data  can  be  found.  The  remaining  differences  are  well  explained 
by  the  different  methods  and  should  not  been  further  interpieted. 


VI.  Conclusions 


Estimation  techniques  to  examine  seismic  network  detection  capabilities  including 
multiwave  detection  criteria  are  well  established.  For  a  representative  global  50-station 
network  -  composed  by  an  international  group  of  seismologists  (CCD/558,  1978)  and 
further  studied  in  this  report  -  magnitude  thresholds  requiring  a  90  per  cent  probability  of 
at  least  four  detecting  stations  range  from: 

3.4^m^3.6  for  Scandinavia  and  Europe 
3.7^ ^4.0  for  North  America,  Asia,  and  Arctica 
3.8^m^4.3  for  South  America,  Africa,  and  Antarctica 
4.0^mj^4.6  for  Australia,  New  Zealand,  and  Pacific. 

Current  telescismic  detection  thresholds  are  at  least  half  a  magnitude  unit  higher. 
The  discrepancy  between  model  results  and  actual  reportings  can  be  explained  by  adjust¬ 
ing  station  data.  These  data  have  more  influence  on  global  detection  thresholds  than  vari¬ 
ation  of  amplitude-distance  curves  within  reasonable  bounds.  After  estimation  of  noise 
statistics  from  long  observation  intervals  predicted  magnitude  detection  thresholds  and 
empirical  values  are  in  good  agreement. 

Considerable  improvements  of  existing  seismological  instrumentation  and  practice  is 
needed  to  achieve  those  values  predicted  from  the  model  network: 

•  To  lower  the  background  noise  of  stations  we  recommend  installation  of 
seismometers  in  boreholes  or  use  of  miniarrays  which  offer  the  additional 
benefit  of  phase  identification. 


•  Station  noise  should  be  permanently  monitored  by  reporting  noise  values  in 
connection  with  conventional  station  readings. 

•  Automated  detector  algorithms  should  be  implemented  at  all  stations. 

•  Data  centres  should  use  all  relevant  information  for  defining  events,  i.e.,  use 
more  than  p-phases  (especially  pkp)  in  association  process. 

•  Amplitude  attenuation  tables  for  all  distances  should  be  standardized  to 
specific  frequency  bands  and  applied  to  amplitude  measurements  from 
identical  instrument  responses. 

With  these  modifications  a  global  50-station  network  is  assumed  to  come  up  to  telese- 
ismic  detection  thresholds  as  predicted  from  our  calculations. 
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